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ABSTRACT 


A  new,  relatively  simple  and  inexpensive 
non-destructive  inspection  technique  for  advanced 
composite  materials  is  proposed  and  its  feasibility  is 
demonstrated.  This  technique  uses  encapsulated 
cholesteric  liquid  crystals  to  sense  small  temperature 
differences  which  result,  from  nonuniform  heat  transfer 
through  composites.  Discontinuities  in  heat  transfer 
evidenced  by  contrasting  surface  color  patterns 
indicate  material  discontinuities,  i.e.  flaws. 

Preliminary  investigations  into  the  thermal 
conductivity  of  a  .04  1  in.  thick  graphite/epoxy 
laminated  panel  in  the  direction  normal  to  the 
composite  laminae  and  in  the  direction  parallel  to  the 
laminae  are  described.  The  coefficients  of 
conductivity  in  the  two  directions  were  found  to 
differ  by  an  order  of  magnitude.  A  qualitative  test 
of  the  technique  in  locating  embedded  teflon  triangles 
is  reported,  and  a  design  for  a  testing  device  is 
proposed. 
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LIST  OF  SYMBOLS 


A  the   cross-sectional   area   of  a  fin  or  wall  through 
which  heat  flow  is  examined  —  ft2 

C  the   circumference   of   a  cross  section  at  a  location 
along  the  length  of  a  fin  —  ft 

H  a  constant  —  h  /km  —  Nondimensional 

e 

h  lateral  fin  surface  film  coefficient  —  BTO/h  ft2  °F 

h  thermal  contact  coefficient  —  BTO/h  ft2  °F 

c 

h  film  coefficient  at  fin  end  —  BTO/h  ft2°F 

e 

k  thermal  conductivity  —  BTO/h  ft  °F 

k  thermal  conductivity  of  material  between  planes  m 

mn  and  n  —  BTO/h  ft  <>F 

L  fin   length   measured   perpendicular   to   the    major 
surface  of  attachment  —  ft 

m  a  constant  —  the  square  root  of  hC/kA  —  ft  * 

g  rate  of  heat  flow  through  a  material  —  BTO/h 

t  temperature  —  °F  or  °C 

t  temperature  of  the  bulk  of  ambient  fluid  —  °F  or  °C 
f 

t  temperature  at  location  n  —  °F  or  °C 
n 

t  fixed  temperatute  at  the  fin  base  —  °F  or  °C 
o 

x  coordinate  along  the  length  of  the  fin  —  ft  or  in 

Ax  material  thickness  —  ft  or  in 

Ax  material  thickness  between  planes  m  and  n  --  ft  or  in 

mn 

d  temperature   difference  between  a  point  in  a  material 
and  the  ambient  fluid  —  °F  or  °C 

9       temperature   difference  between  fin  base  temperature 

o  and  the  ambient  fluid  —  °F  or  °C 

1/b  A  thermal  contact  resistance  —  h  ©F/BTO 
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INTRODUCTION 


Sines  the  early  development  of  boron/apoxy  and 
graphite/epoxy  composite  laminates  it  has  been  recognized 
that  non-destructive  testing  techniques  developed  for  use  in 
metalic  structures  are  often  inadequate  for  composites.  An 
effort  to  obtain  effective  non-destructive  '  testing 
techniques  has  resulted  in  the  development  of  a  number  of 
moderately  successful  methods.  These  techniques  include 
ultrasonic,  radiographic,  sonic,  infrared,  thermochromic, 
photochromic,  holographic  interf erometr y,  dye  penetrant  and 
hardness  tests,  and  are  discussed  in  Refs.  1  and  2. 
Although  each  method  has  its  unique  advantages  and 
disadvantages,  they  generally  share  a  high  level  of 
sophistication  and  complexity  and  frequently  enjoy  only 
limited  laboratory  use.  It  is  not  surprising,  therefore, 
that  a  simple  yet  effective  non-destructive  testing 
technique  is  desired  which  can  be  inexpensively  and  easily 
performed  by  personnel  with  limited  technical  training. 

In  recent  years  some  potentially  effective  ultrasonic 
and  holographic  techniques  have  been  devised  (Ref.  3)  and 
the  use  of  a  relatively  compact  mobile  holographic 
non-destructive  testing  device  has  been  demonstrated 
(Refs.  4  and  5) .  Among  the  simplest  are  those  developed  in 
the  latter  part  of  the  1960's  which  involve  the  use  of 
cholesteric  liquid  crystals  in  thermal  mapping.  Some 
applications  of  these  methods  are  presented  in 
Refs.  6  and  7,  including  the  location  of  interlaminar 
electrical  shorts  in  electronic  circuit  boards,  the 
prediction  of  fatigue  and  fracture  locations,  crack 
detection   and   the   inspection   of    adhesively    bonded 


structures.  Brilliant  color  patterns  can  be  observed  through 
the  application  of  liquid  crystals  to  a  few  square  feet  of  a 
surface  which  exibits  nonuniform  temperatures.  Their  use  is 
particularly  desirable  when  the  temperature  variation  can  be 
attributed  to  nonuniform  heat  transfer  through  a  material 
containing  some  nonhomogeneity  or  discontinuity.  Two 
problems  that  hindered  the  early  use  of  cholesteric  liquid 
crystals  in  large-scale  non-destructive  testing  are  the 
short  useful  lifetime  of  pure  cholesteric  liquid  crystals 
and  the  expense  of  coating  large  objects.  The  following 
development  of  a  non-dsstructive  technique  for  detecting 
discontinuities  and  nonhomogeneities  in  advanced  composite 
plates  addressed  the  problem  of  coating  large  areas. 

Three  separate  stages  in  the  development  of  a  proposed 
compact,  mobile  device  were  completed.  First  was  the 
determination  of  the  coefficient  of  thermal  conductivity 
normal  to  the  layers  of  a  graphite/epoxy  panel.  Second  was 
the  determination  of  the  coefficient  of  conductivity 
parallel  to  layers.  References  8  and  9  provided  the 
theoretical  background  for  a  steady  state  analysis  of  heat 
transfer  in  the  laminate.  When  order-of-magnitude  values  for 
the  two  coefficients  of  conductivity  had  been  obtained,  the 
third  stage  was  begun.  It  involved  the  actual  test  of  a 
14  x  14  in.  graphite/epoxy  panel  containing  intentionally 
embedded  nonhomogeneities.  Ultimately  the  final  device  was 
conceptually  designed  and  partially  constructed.  Key  to  the 
design  was  the  procurement  of  a  thin  (.005  in.) 
mylar-insulated  foil  heater  6  in.  square. 

As  previously  mentioned,  pure  liquid  crystals,  while 
exibiting  brilliant  colors,  deterorate  rapidly  once  applied 
to  a  surface  and  provide  only  a  few  hours  for  viewing  and 
experimentation.  A  brief  discussion  of  liquid  crystals 
including  encapsulated  cholesteric  liquid  crystals  (the  type 
actually   used   in   the   work   hereinafter    described)    is 
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appropriate. 

1 •   Cholesteric  Liquid  Crystals 

Detailed  information  concerning  the  origin, 
composition  and  uses  of  liquid  crystals  can  be  found  in 
Refs.  10,  11  and  12,  and  is  summarized  here.  In  the  latter 
part  of  the  19th  and  early  20th  centuries,  scientists  such 
as  Friedrich  Reinitzer  began  work  on  selected  organic 
compounds  that,  in  the  thermal  transformation  from  solid  to 
liquid  state,  passed  through  an  intermediate  phase  which 
exibited  the  anisotropic  optical  properties  of  crystalline 
solids.  Subsequent  research  concerning  the  structure  of 
liquid  crystals  resulted  in  classification  of  liquid 
crystals  into  one  of  three  categories:  smetic,  nematic,  or 
cholesteric,  the  category  being  determined  by  the  particular 
molecular  arrangement  of  the  liquid  crystals. 

The  response  of  cholesteric  liquid  crystals  to 
thermal  variations  makes  them  an  important  tool  for  the 
acquisition  of  both  quantitative  and  qualitative  information 
in  laminate  heat  transfer  problems.  A  particular  formulation 
of  specific  cholesterol  esters  will  provide  liquid  crystals 
that  exibit  a  reversible  and  repeatable  response  to 
temperature  of  a  given  range.  The  cholesteric  liquid 
crystals  respond  by  sequentially  passing  through  the 
complete  visual  spectrum  —  red  through  violet  --  in 
fractions  of  degrees,  or  multiples  of  degrees,  depending  on 
the  specific  formulation.  The  range  is  generally  referred  to 
as  the  "event  temperature  range"  or  "colorplay  range". 
Since  colors  scattered  by  liquid  crystals  represent  only  a 
fraction  of  the  incident  light  (the  remainder  being 
transmitted  by  the  liquid  crystals) ,  an  absorptive  black 
background  is  needed  to  prevent  reflection  of  the 
transmitted  light. 
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Many  of  the  problems  posed  by  the  use  of  pure 
cholesteric  liquid  crystals,  including  atmospheric 
contamination  and  short  life,  have  been  greatly  reduced  by 
an  encapsulation  process  for  packaging  labile  materials 
developed  by  the  National  Cash  Register  Company.  The 
encapsulation  process  coats  the  liquid  crystals  with  gslatin 
in  a  polyvinyl  alcohol  binder  to  form  tiny  packages  2-50 
microns  in  diameter.  In  addition  to  extending  the  life  of 
the  liquid  crystals  to  as  long  as  several  years  by 
minimizing  surface  contamination  and  ultra  violet  light 
exposure  of  the  raw  liquid  crystals,  the  encapsulation 
process  greatly  reduces  the  variation  in  color  due  to 
viewing  angle.  The  colors  remain  readily  visible,  albeit 
subdued,  to  the  human  eye.  Encapsulated  Liquid  Crystals  may 
be  obtained  in  two  forms,  either  precoated  on  a  blackened 
paper  or  mylar  film,  or  in  a  water  based  slurry.  The  slurry 
form  was  used  for  spray  coating  surfaces  after  thinning,  1 
part  E.L.C.  and  2  parts  water,  with  the  addition  of  2  drops 
of  Kodak  Photoflo  solution  to  2  oz.  of  the  mixture  to 
improve  surface  wetting. 

Encapsulated  Liquid  Crystals  may  be  formulated  to 
respond  over  a  wide  range  of  event  temperatures.  There  are 
three  series  represented  by  E.L.C.  nomenclature: 
5  (sensitive),  1.5  ♦  0.5°C  colorplay;  R  (regular), 
3.5  +  1.0°C  colorplay;  W  (wide),  6.0  +  2.0<>C  colorplay. 
The  S  -  series  were  exclusively  used  '  for  highest 
sensitivity;  their  characteristic  event  temperatures  are 
listed  in  Appendix  A. 
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II.   HEAT  TRANSFER  NORMAL  TO  LAYERS  OP  A  LAMINATE 


ONE-DIMENSIONAL  STEADY  STATE  HEAT  CONDUCTION 


"Steady  state  conduction"  is  a  condition  which  exists 
when  temperatures  of  fixed  points  within  a  heat-conducting 
body  do  not  change  with  time.  This  also  implies  that  the 
time  rate  of  heat  flow  between  any  two  selected  points  is 
constant  with  time.  The  term  "one-dimensional"  can  be 
applied  to  a  heat  conduction  problem  when  only  one  space 
coordinate  is  necessary  to  describe  the  temperature 
distribution  within  a  heat-conducting  body.  Although  such  a 
situation  rarely  exists  in  practice,  the  one-dimensional 
condition  can  be  assumed  as  a  close  approximation  in  the 
study  of  heat  flow  normal  to  the  surface  of  the 
graphite/epoxy  laminate.  The  heat  flow  through  a  plane 
wall,  at  regions  sufficiently  removed  from  the  edges, 
depends  only  on  the  coordinate  measured  normal  to  the  plane 
of  the  wall.  By  insulating  a  finite  wall  at  its  edges  one 
may  neglect  edge  effects  and  a  practical  one-dimensional 
heat  conduction  problem  exists. 

The  elementary  heat  conduction  problem  of  a   plane   wall 

of   finite   thickness  Ax  but  infinite  in  extent  (  to  allow  a 

one-dimensional  assumption)  with  each  face  maintained   at   a 

uniform  temperature  is  represented  by  Eg.  2.1  (Ref.  8)  . 

g      t  -  t       t  -  t 

_L  =  k  _i 2_=  jc  i L.  (2.1) 

A      x  -  x        Ax 

2     1  12 

The  above  expression  is  written  in  terms  of  g/A,  the  rate  of 
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heat  flow  per  unit  area,  Ax,  the   wall   thickness,   and  the 
respective  uniform   face  temperatures,   t   and   t  ,  for  a 

12 

material  having  thermal  conductivity  k. 

Mow  consider  the  more  general  plane  wall,  illustrated  in 
Fig.  1  (a) ,  composed  of  three  layers  (although  any  number 
could  have  been  used)  each  of  different  thickness  and 
thermal  conductivity.  Denote  each  juncture  of  two  different 
materials  by  a  number  1  through  4  and  adopt  the  following 
notation : 


■Ax    =  thickness  between  planes  m  and  n 
mn 

■   k    =  thermal  conductivity  of  material  between 
mn 

planes  m  and  n 


If  the  outside  temperatures  ,  t   and  t  ,   are   specified 

and    the    individual    material  thicknesses  and   thermal 

conductivities  are  known,  then  the  temperature  distribution 
through  the  wall  and  rate  of  heat  flow  can  be  calculated. 
Since  the  steady  state  condition  is  assumed  to  exist  and 
since  the  rate  of  heat  flow  per  unit  area  q/A,  is  constant, 
Eg.  (2.1)  may  be  applied  to  each  layer  yielding  the 
following: 

q  t   -    t  t-t  t-t 


1  2         =  2  3         =  3 


Ax      /k  Ax      /k  Ax      /k 

12  12  23         23  34  34 


(2.2) 


Solving   for  the  individual  temperature  differences  and  then 
eliminating  the  unknown   temperatures   the   expression,   Eg. 
(2.3),   for   heat   flow  through  the  wall  in  terms  of  overall 
temperature  difference  is  obtained. 
5         t  -  t 


A    Ax     Ax     Ax 

1  2   +     23   +     34 

k      k      k 

12        23        3  4 


(2.3) 


Next  consider  a  practical  modification  to  the   "infinite 
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slab"  (the  wall  with  no  edge  effects)  problem  just 
discussed.  By  extracting  a  small  circular  section  normal  to 
the  plane  of  the  wall  in  Pig.  1 (a)  and  insulating  it,  the 
one-dimensional  problem  can  be  preserved,  with  the 
assumption  that  although  edge  effects  may  now  exist,  they 
are  small.  The  practical  problem  now  depicted  in  Fig.  1(b) 
is  one  of  two  insulated  circular  cylinders  of  the  same 
material  sandwiching  a  thinner  dissimilar  disc. 

1 .   Contact  Resistance 

Imagine  two  solid  rods  brought  into  contact  as 
indicated  in  Fig.  2(a),  with  sides  insulated  so  that  heat 
flows  only  in  the  axial  direction.  The  materials  may  have 
different  thermal  conductivities,  but  if  their  sides  are 
insulated  the  heat  flux  must  be  the  same  through  both 
materials  under  the  steady  state  conditions.  Experience  has 
shown  that  the  actual  temperature  profile  through  the 
materials  varies  approximately  as  shown  in  Fig.  2(b).  The 
temperature  drop  at  plane  2  is  caused  by  thermal  contact 
resistance.  The  physical  mechanism  of  contact  resistance  may 
be  better  understood  by  close  examination  of  the  joint  in 
detail.  No  real  surface  is  perfectly  smooth  and  the  actual 
surface  roughness  is  believed  to  play  a  central  role  in 
determining  the  contact  resistance.  There  are  two  principal 
contributions  to  the  heat  transfer  at  the  joint: 

■  The  solid-to-solid  conduction  at   the   spots   of 
contact 

■  The  conduction  through  entrapped   gases   in   the 
void  spaces  created  by  the  contact. 

The   second   factor   represents  the  major  resistance  to  heat 

flow.  The  quantity  1 / (h  A)  is  the  thermal  contact  resistance 

c 

and   h    is   the  contact  coefficient.   Including  the  contact 
c 
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coefficient  in  Eq.  (2.3)  accounts  for  the   thermal   contact 

resistance   encountered  in   the  real  problem  and  Eg.  (2.4) 
results. 

q  t  -  t 

=  1   _4 


A   Ax     Ax     Ax     _L    ,  1  (2.4) 

k      k      k      c     c 

12        23         34        23       34 

Let  the  materials  between  planes  1  and  2  and   planes 
3  and  4  be  the  same  and  their  thicknesses   x     and   x   be 

12  34 

equal.  Further,  if  the  surfaces  of  contact  at  planes  2  and 
3  have  equal  roughness  characteristics,  one  may  combine  the 
contact  coefficients  and  Eq.  (2.4)  can  be  simplified  as 
follows: 

a  t :  -  + 

(2.5) 


where 


g 

t     -     t 
1               4 

A 

2Ax 

1  2 

k 

12 

Ax             _JL 

♦           23     +  TT~ 

k               c 

23                     2  4 

T^-  = 


h     +  h  (2.6) 

c     c      c 

24       2  3         34 

Once  again  considering  the  experimental  problem  of 
the  thin  disc  sandwiched  between  two  circular  rods 
previously  discussed  and  illustrated  in  Fig.  1(b),  q/A  may 
be  found  using  either  the  first  or  third  of  Eqs.  (2.2). 
Further,  given  the  physical  dimensions  of  the  two  rods  and 
the  disc  as  well  as  their  respective  thermal  conductivities, 
the  contact  coefficient  may  be  found  by  solving  Eg.  (2.5). 

Assume  that  the  contact  conditions  for  two  discs, 
one  of  known  thermal  conductivity  and  one  of  unknown  thermal 
conductivity,  can  be  made  equal.  The  contact  coefficient 
found  in  the  known  case  may  be  used  in  the  other  case  to 
solve  for  the  unknown  thermal  conductivity.  The  problem 
studied  was  one  in  which  the  thermal  conductivity  of  the 
graphite/epoxy   laminate  disc  was  unknown  in  addition  to  the 
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unknown  contact  coefficient. 
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(a)    One-dimensional,   multilayer     wall      case 


(b)  One-dimensional,  thin   specimen    case 


Figure  1  -   HEAT  TRANSFER  NORMAL  TO  THE  LAYERS  OF  A  LAMINATE 
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(a)  Two    insulated    solid   rods 


(b)  Temperature     Profile 


Figure  2  -    CONTACT  RESISTANCE 
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B.   EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  apparatus  used  foe  the  experimental  determination  of 
the  coefficient  of  conductivity  normal  to  the  layers  of  a 
graphite/epoxy  specimen  is  shown  in  Figs.  3  and  4.  Heat 
Transfer  Unit  I,  shown  in  detail  in  Pig.  3,  is  composed  of 
three  basic  components: 

1.  Bottom  boiler  —  steam  heat  source 

2.  Connecting  insulated  copper  rods 

3.  Top  insulated  ice  bath  --  heat  sink 

The  test  section  to  hold  specimen  discs  is  located  in  the 
gap  formed  between  the  ends  of  the  two  rods. 

Each  copper  rod  was  fitted  with  four  small 
(approximately  .020  in.  diameter)  copper-constantan 
thermocouples  inserted  into  holes  drilled  to  the  rod  center 
at  locations  consecutively  rotated  90°. 

The  entire  apparatus,  shown  in  Fig.  4,  consists  of  four 
main  components: 

1.  Voltmeter   capable   of   millivolt   output    for 
tnermccouple  reading 

2.  3-position  rotary  selector  switch 

3.  Ice  bath  reference  point 

4.  Heat  Transfer  Unit  I 

Appendix  B  contains  information  regarding  the  dimensions  and 
thermal  conductivities  of  material  specimens  used,  in 
addition  to  properties  of  other  materials  used  in  the 
construction  of  the  apparatus. 

The  procedure  described  below  was  followed  to  obtain  a 
practical  estimate  of  the  contact  coefficient  which  was  used 
in  Eq.  (2.5)  and  ultimately  in  conjunction  with  Eg.  (2.2)  to 
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find  the  thermal  conductivity  normal  to  the  layers  of  the 
graphite/epoxy  composite.  During  one  testing  period  nine 
3/4  in.  diameter  specimen  discs,  six  metal  and  three 
composite,  were  examined.  Also  during  each  period  two  test 
sequences  were  performed  with  no  specimen  to  determine  the 
magnitude  of  the  copper-heatsink  compound-copper  contact 
coefficient.  For  all  contact  surfaces,  already  with 
similarly  machine  finished  surfaces,  G  and  c  Electronics 
heatsink  compound  was  used  as  a  thin  uniform  coating  to 
reduce  the  effect  of  air  entrapped  in  surface  voids. 

Two  methods  were  used  in  the  arrangement  of  thermocouple 
output: 

•  "The  direct  method,"  which  included  an  ice  bath 
reference  in  the  thermocouple  circuit. 

■  "The  differential  method,"  which  eliminated  the 
ice  bath  by  connecting  the  copper  sides  of  pairs 
of  thermocouples  located  equidistant  from  the 
gap  (1  and  8,  2  and  7,  3  and  6,  4  and  5)  . 

The  following  is  a  detailed  sequence  describing  one  test 
period: 

1.  After  lightly  coating  the  contact  surfaces  of 
one  specimen  disc  with  the  heatsink  compound  it 
was  inserted  in  the  gap  between  the  copper  rods 
with  boiling  water  in  the  bottom  container  and 
the  ice  bath  in  the  top. 

2.  A  25  lb.  weight  was  placed  on  top  of  the  ice 
bath  to  enhance  the  release  of  entrapped  air  in 
the  contact  voids. 

3.  When  the  voltmeter  output  settled  indicating 
that  temperatures  in  the  copper  rods  had 
stabilized  at  steady  state  values, the  millivolt 
output  was  recorded  (using  either  the  pre-wired 
direct  or  differential  method)  using  the   rotary 
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switch    to    sequentially    select   specific 
thermocouple  circuits. 

4.  The  temperatures  (or  temperature  differences) 
were  plotted  versus  rod  length.  A  linear 
polynomial  fit  and  extrapolation  were  used  to 
determine  t   and  t   (or  At)  across  the  gap.  (See 

Fig.  5  as  an  example) . 

5.  Steps  3  and  4  were  repeated  three  times  for 
each  specimen  during  each  test  sequence. 

6.  Steps  1  through  5  were  repeated  for  all  nine 
specimens  and  the  two  previously  mentioned  "no 
specimen"  sequences  were  performed  to  complete 
an  entire  test  period. 

Subsequent  data  reduction  for  appropriate  temperatures 
and  use  of  the  known  thermal  conductivities  provided 
solutions   to   Eq.   (2.5)   in  terms  of  h  .  The  mean  of  these 

contact  coefficients  was  then  used  to   compute   the   thermal 

conductivity   in   a   direction   normal   to  the  layers  of  the 
graphite/epoxy  specimens. 
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Figure  3  -   HEAT  TRANSFER  UNIT  I 
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Figure    5    -       EXAMPLE   At    FOR    TWO    SPECIMENS 
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RESULTS 


Figure  6  is  a  summary  of  the   values   of   1/h    obtained 

c 

over   three   separate   test   periods   displayed   for  viewing 

convenience  and  not  to  infer  a  correlation  between   specimen 

thickness  and  contact  coefficient. 

One  of  the  three  graphite/epoxy  specimens  had   smoother 

surfaces   than   the  other  two;  therefore,  in  calculating  the 

coefficient  of  conductivity   using   data   from   the   rougher 

specimens,   the   value   of   1/h   used  was  increased  from  the 

c 

mean  by  one  sample  standard  deviation   to   account   for   the 

effect  of  increased  surface  roughness. 

The  nine  subsequent  values  for  the  coefficient  of 
thermal  conductivity  measured  normal  to  layers  in  the 
graphite/epoxy  laminate  ranged  from  .534  to. 907  BTU/h  ft  °F 
with  a  sample  mean  value  of  .748  BTU/h  ft  °F  and  associated 
sample  standard  deviation  of  .121  BTU/h  ft  °F  . 
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III.   HEAT  TRANSFER  PARALLEL  TO  LAYERS 


A.   STRAIGHT  FIN  OF  UNIFORM  THICKNESS 

When  examining  the  inplane  heat  transfer  problem  one  is 
led  to  tne  study  of  heat  transfer  for  a  thin  fin.  If  the  fin 
proportions  used  in  practice  are  such  that  the  length  of  the 
fin  (the  dimension  measured  normal  to  the  primary  surface  of 
attachment)  is  large  compared  with  its  maximum  thickness, 
one  may  assume  that  the  temperature  of  the  fin  depends  only 
on  a  single  coordinate  measured  in  the  direction  of  the 
fin's  length.  Thus,  when  the  above  simplification  can  be 
made,  the  problem  for  the  solution  of  the  temperature 
distribution  and  heat  flux  rate  becomes  a  one-iimensional 
problem.  The  following  analysis  addresses  only  the  steady 
state  case. 

For  the  straight  fin  of  uniform  thickness  the 
cross-sectional  area  for  heat  flow  (normal  to  the  length 
coordinate)  is  constant.  Also  the  exposed  surface  area  for 
heat  convection  is  a  linear  function  of  the  distance 
measured  along  the  length;  i.e.,  the  perimeter  of  the  cross 
section  is  constant. 

Figure  7  (a)  illustrates  the  fin  configuration  and  shows 
the  notation  to  be  used  in  the  following  analysis.  The 
symbol  L  denotes  the  length  of  the  fin,  k  the  thermal 
conductivity,  h  the  film  coefficient  of  the  exposed  surface, 
t   tne  fixed   temperature   at   the   fin   base,   and   t   the 

3  f 

temperature    of    the   bulk   of   the   ambient   fluid.    The 

coordinate  distance  along  the  fin  length  is  symbolized  by  x. 
A   and   C   are   used   to   denote  the  area  of  a  cross  section 
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normal  to  x  and  the  perimeter  of  the  section,  respectively. 
From  conservation  of  energy  in  the  steady  state, 
Fourier's  law  and  the  definition  of  the  film  coefficient  h, 
one  obtains  the  classic  differential  equation  for  the 
distribution  of  temperature  in  the  fin  (Ref .  8) : 


dzt    hC 

(t  -  t  )  =  0 

dx2    k\  f 


(3.1) 


Equation  (3.  1)  must  be  satisfied  at  every  point  in  the  fin, 
and  if  the  film  coefficient  and  the  thermal  conductivity  may 
be  assumed  constant,  it  can  be  easily  solved.  The  equation 
can  be  made  more  concise  by  defining  the  following 
parameters: 


0  =  t 


(3.2) 


the 


(3.3) 

differential 


0 


(3.4) 


Substituting   the   definitions    yields 
equation,  Sq.   (3.4). 

d28 

dx2 

The   general   solution   to   Eg.   (3.4)   including   arbitrary 

constants  B  and  D  found  by  applying  boundary   conditions   is 

given  by  Eq .  (3.5). 

-mx     -mx 
9  =  Be    +  De  (3.5) 

For   the   present   1st   the  boundary  conditions  be  those 

depicted  in  Fig.  7(a)  --  in  imposed  temperature  at  the   base 

and   convection   at   the  free  end.  These  boundary  conditions 

may  be  written: 

e  =  e  =  t  -  t 

0     0     f 

(3.6) 


At  x  =  0 
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At  x  =  L:  = 


h  9 
dx     Tc 


It  should  be  noted  that  the  boundary  condition  of  convection 

at   the   end   of   the   fin  as  stated  above  permits  the  film 

coefficient  at  the  end,  h  ,  to  be  different-  from   the   film 

e 
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coefficient    on    the    lateral   surfaces,   h.    Subsequent 

evaluation  of  B  and  D   and   substitution   into   the   general 

solution    leads    to    the   following   expression   for   the 

temperature  distribution  along  the  fin: 

9    t  -  t    cosh  m  (L-x)  +  Hsinh  m  (L-x) 

t-  ~ T— 7 ~  (3-7) 


9    t  -  t        cosh  mL  +  Hsinh 
o    o    f 


where 


h  /hC 

H  =  _£.    and     m  =-i/ (3.8) 

km  ykA 

In  the  development  of  Eq.  (3.7)  it  was  assumed  that  the 

fin   was   slender   enough  for  a  one-dimensional  condition  to 

prevail.  If  this  conditioa  is  met  it  is  very  likely  that  the 

heat  convected  out  the  end  of  the  fin  is  a  small  fraction  of 

the  heat  convected  out  the  other   surfaces.   Therefore,   Eq. 

(3.7)   may   be   simplified  by  neglecting  end  heat  loss.  This 

implies  h   =0,  or  H=0  and  Eq.  (3.7)  becomes: 
e 

9  t  -  t     cosh  m  (L-x) 

-  ■  L  =    ——  (3.9) 

9  t  -  t       cosh  mL 

o  o    f 

A  major  difficulty  in  using  Eg.   (3.9)   in   a   practical 

problem   is   in   the   determination   of   h.  The  basis  for  an 

experimental  determination  of  h  is  illustrated  in  Fig.  7(b), 

in   which   the   temperature  distribution  along  the  length  of 

two  fins  is  known  as  well   as   their   dimensions.   If   these 

temperature   distributions  were  derived  from  the  same  actual 

conditions  for  two  fins  of  identical  dimensions,   both   fins 

must   have   the  same  film  coefficient  h.  By  suostituting  the 

appropriate  known  values  into  Eq .  (3.9)  for  each  of  the   two 

fins   and   ratioing   the   subsequent   equations,   Eg.  (3.10) 

results. 

9      9  (t-t    )     (t    -t    )  cosh    m    (L-x)     cosh    m    L 

_1 02      =  f      1  0         f     2      =     X I  2      2  (3.10) 

9      9  (t-t    )     (t    -t    )  cosh   m    L      cosh    m    (L-x) 

012  f     2  0         p      l  l      l  2  2 

Since  the  only  unknown  is  now  h,  common   to   both   fins   and 
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found   in   the   terms   m   and  m  .  an  iterative  solution  will 

1  2 

yield  the  actual  value  of  h.  Once  the  actual  value  of   h   is 

known  it  may  be  used  to  solve  for  an  unknown  coefficient  of 
thermal  conductivity  in  any  other  geometrically  similar  fin 
sharing  the  same  environment. 
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(a)    Straight    fin   of  uniform    thickness 


(b)  Temperature    vs.  Length   —  two    fins 

Figure  7  -   STRAIGHT  FINS 
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B.   EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

To  obtain  information  concerning  the  thermal 
conductivity  parallel  to  layers  in  samples  of  specific 
graphite/epoxy  laminates,  Heat  Transfer  Unit  II  was 
constructed  (Fig.  8).  It  consists  of  a  rectangular  aluminum 
tank  with  thres  tubular  threaded  portions  welded  into  holes 
in  the  tank  side  wall.  Each  of  the  tubular  portions  was 
fitted  with  an  open  center  cap  and  teflon  washer.  Six  fin 
units  were  used;  each  was  a  long,  narrow,  thin  fin,  inserted 
into  a  size  no.  9  rubber  stopper.  For  each  test  period 
three  of  these  fin  units,  two  with  known  thermal 
conductivities  and  one  with  unknown  thermal  conductivity, 
were  placed  in  the  tubular  sections  and  caps  tightened  to 
lock  the  units. 

All  of  the  fins  were  aachined  as  nearly  as  possible  to 
the  same  dimensions,  and  all  were  prepared  for  testing  in 
the  same  manner.  A  typical  fin  unit  is  shown  in  Fig.  9.  The 
base  of  each  fin,  (considered  to  be  the  distance  from  the 
exterior  face  of  the  stopper,  through  the  stopper  to  the 
interior  of  the  tubular  section  of  the  boiler  tank) ,  a 
distance  of  3  in.,  was  coated  with  a  clear  plastic 
varathane.  This  coating  was  primarily  necessary  for  the 
waterproofing  of  the  laminate  fins.  The  portion  of  the  fin 
exterior  to  the  stopper  over  which  a  temperature  profile  was 
to  be  measured  was  sprayed  with  flat  black  enamel.  Over  the 
enamel  was  sprayed  a  mixture  of  cholesteric  liquid  crystals, 
S-30,  S-32,  and  S-62. 

During  a  single  testing  period  water  in  the  rectangular 
tank  was  brought  to  a  boil  with  the  use  of  an  alectric  hot 
plate  which  ultimately  provided  a  steady  state  temperature 
base  for  the  fins  in  well  circulated  boiling  water. 

Whan  the  color  play  patterns  on  the  three  fins  had 
stabilized,  the  ambient  temperature   was   rscorded   and  the 
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positions  of  the  green  and  blue  color  zones  were  marked  on 
each  fin  as  well  as  the  position  of  the  black  separation 
zone  between  s-30  and  s-32  crystals.  The  position  of  each 
mark  was  measured  from  the  face  of  the  stopper  and  recorded 
with  the  appropriate  temperature  for  the  particular  color 
zone.  The  length  of  the  fin  extending  from  the  face  of  the 
stopper  was  also  recorded. 

From  these  measurements  temperature  vs.  fin  length 
curves  were  constructed.  Figure  10  is  a  composite  plot  of 
six  individual  curves  found  in  Appendix  C.  Incrementing  x 
from  0.5  in.  to  2.0  in.  in  0.5  in.  steps  and  using 
corresponding  temperatures  from  fins  of  known  thermal 
conductivity,  a  value  for  h  was  determined.  This  value  of  h 
was  used  to  calculate  a  value  for  k  measured  parallel  to 
layers  in  the  two  laminate  fins.  A  typical  representation  of 
these  values  of  thermal  conductivity  is  given  in  Fig.  11. 
The  value  for  the  thermal  conductivity  of  the  graphite/epoxy 
laminate  in  a  direction  parallel  to  layers  was  taken  to  be 
the  average  of  the  values  represented  in  Fig. 11  for  each 
test  period. 
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C.   RESULTS 

A  summary  of  the  resulting  values  for  the  coefficient  of 

conductivity   in   a   direction   parallel   to   layers   in  the 

graphite/epoxy  fins  is  presented  in  Fig.  11.   The   scattered 

results   at   1.5  and  2.0  inches  are  attributed  to  t  -  t  the 

f 

sensitivity  of  the  calculations  to  small  values  of  (Fig.  10) 

obtained  in  that  region.  Including  the  scattered  values, 
however,  still  yields  a  sample  mean  value  for  k  parallel  to 
layers  in  the  graphite/epoxy  fin  of  6.6  1  BTU/h  ft  °F  with 
sample  mean  value  of  1.25  3TU/h  ft  °F. 
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IV.   COMPOSITE  TEST  PANEL 


A.   BACKGROUND 

With  knowledge  of  the  relative  magnitudes  of  thermal 
conductivities  for  the  specific  gr aphite/epoxy  laminates 
available,  the  final  step  before  designing  an  actual  flaw 
detection  device  was  to  attempt  flaw  detection  in  a 
graphite/epoxy  panel  with  intentionally  embedded  flaws. 

Recalling  the  discussion  of  one-dimensional  heat  flow 
through  a  wall,  and  observing  that  the  panel  could  be 
treated  as  a  thin  wall,  provides  a  basis  for  a 
non-destructive  test  of  the  panel.  Consider  two  paths  for 
heat  flow  from  the  rear  surface  to  the  front  surface  of  the 
panel.  Let  Path  1  pass  through  only  the  wall  material  and 
Path  2  pass  through  the  wall  material  as  well  as  through  an 
embedded  flaw.  3y  supplying  a  uniform  rear  surface 
temperature  and  assuming  heat  flow  rate  constant,  and  with 
the  application  of  Eg.  (2.3)  to  both  paths,  it  can  easily  be 
shown  that  the  front  surface  temperature  will  be  different 
for  the  two  paths.  Is  the  temperature  difference  between  a 
surface  location  above  a  flaw  a«nd  a  surface  location  below 
which  no  flaw  exists  sufficient  to  be  detected  through  the 
use  of  encapsulated  cholesteric  liquid  crystals?  That  is 
the  Question  addressed  in  the  experiment  described  below. 


3.   APPARATUS  AND  PROCEDURE 

An   3-ply   [ (0/90)  ]    graphite/epoxy   panel,   shown   in 
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Fig.  12,   was  constucted   from   NARMCO   Rigidita    5208-QCC 

Thornel  300  prepreg  tape.  Sixteen  teflon  triangles,  shown  in 
detail  in  Fig.  13,  were  embedded  in  the  panel  as  intentional 
"flaws"  or  material  discontinuities. 

Three  types  of  quantitative  information  regarding  the 
detection  sensitivity  were  desired: 

■  Flaw  width  sensitivity  information 

■  Flaw  thickness  sensitivity  information 

■  Flaw  depth  location  sensitivity  information 

The  area  of  the  triangles  and  the  included  angles  were 
selected  to  facilitate  evaluation  of  the  anticipated  surface 
detection  pattern.  Choosing  the  three  angles  such  that 
angles  2  and  3  were  twice  and  four  times  angle  1, 
respectively,  allowed  consideration  of  three  substantially 
different  tapers. 

It  was  anticipated  that  by  comparing  the  shape  of  the 
surface  temperature  pattern  with  the  known  shape  of  the 
teflon  triangles  one  could  quantitatively  determine  flaw 
width  sensitivity  of  the  proposed  technique.  The  1  sq .  in. 
triangle  area  chosen  was  simply  to  facilitate  data 
reduction. 

To  obtain  thickness  sensitivity  information  the  teflon 
triangles  were  constructed  with  four  different  thicknesses, 
each  thickness  consecutively  increased  by  the  thickness  of 
the  first.  The  resulting  flaw  thickness  range  extended  from 
6%    to  21^  of  the  finished  panel  thickness. 

As  the  layup  process  proceeded  four  teflon  triangles, 
varying  in  thickness  only,  were  positioned  in  each  of  four 
interlaminar  interfaces,  as  depicted  in  Fig.  12.  By 
observing  the  panel  from  both  sides  one  obtained  data  for 
every  possible  flaw  depth. 

In  order  to  accurately  locate  the  triangles  after  the 
panel  had  been  cured,  the  triangles  were  laid  on  a  carefully 
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measured  grid  which  was  marked  daring  layup  by  embedding 
0.002  in.  wires  in  the  0.5  in.  border  of  the  panel. 

After  curing,  the  panel  surfaces  were  sprayed  with  flat 
black  enamel.  When  the  enamel  had  dried,  the  surfaces  were 
then  sprayed  with  encapsulated  cholesteric  liquid  crystal 
S-30. 

The  problem  of  uniformly  heating  the  panel  was 
approached  by  mounting  the  panel  in  the  vertical  position  in 
a  small  vice  and  attaching  the  mouth  of  a  large  tubular 
polyurethane  plastic  bag  to  the  edges  of  the  panel,  Fig.  14. 
The  attachment  was  achieved  with  long  strips  of  U-shaped 
plastic  molding  which  left  1.5  in.  openings  at  the  corners 
of  the  panel.  In  the  other  end  of  the  bag  was  attached  a  1.5 
in.  diameter  hose.  A  smaller  perforated  plastic  bag  had 
previously  been  attached  to  the  end  of  the  hose  so  that  warm 
air  passing  into  the  larger  bag  would  be  effectively 
diffused.  The  other  end  of  the  hose  was  connected  to  a  small 
heater-blower  unit  capable  of  heating  low  velocity  air  to  as 
high  as  approximately  70°C. 

When  the  heater-blower  unit  was  switched  on,  the  large 
bag  was  slightly  pressurized  with  warm  air  which  soon 
reached  a  near  steady  stata  condition,  thus  providing  nearly 
uniform  heat  over  the  rear  surface  of  the  panel.  3y 
observing  the  front  exposed  surface,  one  could  detect  the 
color  patterns  caused  by  non  uniform  heat  conduction  through 
the  panel  influenced  by  the  embedded  teflon  triangles. 
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(a)  Front    view 

—  6  layer  s    over 

' 

* 

, 

1 — 2  lay  ers  over 

(b)  Edge  view 


Figure  12  -    GRAPHITE/EPOXY  LAMINATE  TEST  PANEL 
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L\=  26° 

Z2=5I° 
L  3  =  103° 

Area  =  1.0    in* 


(Shown  actual  size) 


Type 

Thickness  ,  in 

%d 

A 

.009 

21.4 

B 

.0065 

15.4 

C 

.0045 

10.7 

D 

.0025 

6.0 

*  d  =  finis hed  thickness   of  pane 


Figure  13  -   TEFLON  TRIANGLES 
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C.   RESULTS 

With  the  heating  arrangement  used  it  became  apparent 
that,  holding  surface  temperatures  within  a  2<>C  tolerance  at 
approximately  30<>c  for  the  color  play  zone  of  the  liquid 
crystals  was  not  possible.  However,  under  the  slow 
transient  condition  that  was  actually  established,  fifteen 
of  the  triangles  were  observed  in  their  undistorted  shapes. 
Since  actual  flaw  detection  capability  was  the  primary  goal, 
further  attempts  to  maintain  the  2°C  tolerance  on  the 
surface  temperature  were  not  undertaken. 

To  enhance  the  location  of  the  thinnest  teflon  triangles 
embedded  under  the  most,  layers,  it  was  necessary  to  overcome 
the  effects  of  the  different  thermal  conductivities  normal 
to  each  other  in  the  panel,  which  tended  to  diminish  the 
surface  color  patterns.  Enhancement  of  the  pattern  was 
accomplished  by  pulsed  heating  (or  cooling)  of  the  local 
surface  above  the  flaw  with  an  additional  heater-blower. 

Without  sophisicated  photographic  and  printing 
techniques  it  was  impossible  to  present  black  and  white 
photographic  results  that  would  accurately  represent  the 
subdued  surface  color  patterns  that  were  actually  observed. 
However,  the  blue  and  green  contrasting  pattern  most  easily 
observed  indicated  a  degree  of  sensitivity  sufficient  to 
detect  the  surface  temperature  difference  between  areas 
located  above  the  flawed  and  unflawed'  regions  within 
approximately  0.5°C. 
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V .   CONCLUSIONS 

It  is  significant  to  note  that  with  data  from  a  limited 
number  of  test  periods  it  was  possible  to  effectively  and 
consistently  obtain  order  -  of  -  magnitude  values  for  the 
coefficients  of  conductivity  both  normal  and  parallel  to 
layers  in  a  gr aphite/epoxy  laminate.  It  is  also  important 
to  note  that  without  the  use  of  liquid  crystals  and  with  the 
possible  substitution  of  thermocouples  for  temperature 
measurement,  the  procedure  for  determining  the  coefficient 
of  conductivity  parallel  to  layers  described  in  part  III 
would  be  at  best  extremely  difficult.  Finally,  embedded 
triangular  flaws  ranging  in  thickness  from  approximately  6  % 
to  21  %  of  the  total  laminate  thickness  and  at  every 
interlaminar  depth  from  one  to  seven  were  quite  easily 
detected  in  undistorted,  clearly  observable  shapes. 

In  addition  to  inspection  of  the  panel  with 
intentionally  embedded  flaws,  the  pulsed  heating  technique 
was  successfully  employed  in  another  research  study  to 
obtain  quantative  information  about  impact  damage  on  the 
same  material.  Long  before  normally  observable  damage  had 
occurred  to  the  graphite/epoxy  panels,  the  pulsed  heating 
technique  provided  both  location  and  magnitude  of  the 
growing  damage  zone  as  sequentially  increased  loads  impacted 
the  panel.  Only  the  existence  of  damage  was  studied, 
however,  and  no  attempt  was  made  by  this  investigator  to 
determine  the  exact  nature  of  the  discontinuities,  appearing 
as  a  "bruise"  early  in  the  impact  sequence. 
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VI.   RECOMMENDATIONS 


To  add  to  the  library  of  information  associated  with  the 
family  of  advanced  composites  it  would  be  beneficial  to 
fully  investigate  the  transient  heat  transfer  phenomena  of 
the  multilayered  wall  in  which  the  coefficient  of 
conductivity  varies  with  coordinate  direction  —  the 
composite  laminate  panel. 

The  potential  use  of  encapsulated  liquid  crystals  as  an 
investigative  tool,  particularly  in  the  inspection  of  modern 
composite  aircraft  components,  is  great.  A  particularly 
beneficial  area  for  investigation  would  be  a  correlation  of 
failure  prediction,  as  a  result  of  projectile  penetration 
damage,  with  theoretically  calculable  failure  orginating 
from  machined  holes.  An  investigation  of  this  nature  could 
ultimately  provide  the  ability  to  make  practical  on  site 
damage  assessment  and  repair  inspection  of  aircraft. 

Of  primary  interest  was  the  conceptual  design  of  a 
device  chat  used  the  concepts  investigated  in  part  IV.  Key 
to  the  design  was  the  manufacture  of  a  very  thin, 
transparent,  heating  sheet  that  would  provide  a  warm  contact 
"foot"  on  a  slightly  pressurized  clear  plexiglass  box.  A 
heater  of  this  type  was  designed  and  is  under  construction 
by  MINCO  Products,  Inc.,  Minneapolis,  Minnesota. 
Specifications  for  the  0.005  in.  mylar  and  foil  heater  are 
presented  in  Fig.  15.  It  is  recommended  that  investigation 
of  the  non-destructive  testing  device  pictured  in  Fig.  16 
utilizing  the  sheet  heating  element  be  continued.  The 
proposed  device  includes  a  sprayed  on  layer  of  liquid 
crystals  and  layer  of  black  paint  on  the  outside  of  the 
heater  sheet,  eliminating  the  need  to  coat  the  damaged 
surface.  With  airflow  through  the  box,  interface  contact 
could   be   achieved   between  the  surface  to  be  inspected  and 
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the  inflated  "foot".  When  the  heater  was  switched  on  the 
surface  heat  would  be  controlled  with  a  rheostat  and  while 
viewing  through  the  clear  box  one  would  pulse  cold  air  on 
the  surface,  causing  rapid  transitions  through  the 
characteristic  event  temperature  range.  Personnel  who  have 
only  minimal  training  could  perform  the  inspection  using 
such  a  device  and  obtain  accurate  information  on  damage  to 
composite  aircraft  components,  inexpensively. 


49 


border  for 
connection 


coat 
one  side 
with  E  LC 
and    black 


wire   attachment   tabs 


Thickness*        finished    .005" 

Materiah  etched  foil    heater 

clear    mylar   film     insulator 

Lines:  .03  3"  wide  ,  .03  3"    spacing 

Energy    requirements1 

deliver   5  watt  s/in2 

from    15  6  v,    135  & 

Figure   15  -   DETECTION  DEVICE  HEATING  ELEMENT 
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APPENDIX  A 


CHOLESTERIC  LIQUID  CRYSTALS 


liquid 
crystal 

color  transition  zones  °C 
red     green    blue 

s-30 

30.5 

31.0 

31.6 

s-32 

32.5 

33.1 

33.7 

s-34 

34.7 

35.  1 

35.7 

s-36 

36.4 

36.8 

37.4 

s-38 

38.6 

39.1 

39.7 

s-40 

40,6 

41.  1 

41.7 

s-43 

43.2 

43.9 

44.6 

s-45 

45.1 

45.6 

46.2 

s-62 

62.2 

62.6 

63.1 

(Sujunarizad  from  Refs.  13  and  14) 
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APPENDIX  B 


MATERIAL  PROPERTIES 


Specimen    discs 


Mater  i-al 

k,  3TU/h  ft  oF 

thickness  in 

301  stainless 

9.4 

.085 

Monel 

11.3 

.050 

4130  steel 

27.0 

.063 

7075-T6  aluminun 

70.0 

.049 

Brass 

67.7 

.041 

17-7PH  stainless 

9.3 

.033 

Graphite/epoxy  1 



.038 

Graphite/epoxy  2 



.038 

Graphite/epoxy  3 



.033 

Construction   materials 


Cork 

.024 

Copper 

230 

Foam  rubber 

.01-. 02 

Hard  rubber 

.  1  11-. 13 

Phenolic 

.070-. 090 

Teflon 

.150-. 175 

(Summarized    from    Refs.     13    and    14) 
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APPENDIX  C 


FIN  TEMPERATURE  DISTRIBUTION 


The  fins  that  ware  tested  and  subsequently  listed  in  this 
appendix  were  all  .5  in.  wide  and  .041  in.  thick.   Their 
lengths  varied  as  follows: 

3rass 7.27  in. 

4130  steel 7.20  in. 

17-7PH  stainless. .7. 27  in. 

2024-T3  aluminum. .7. 22  in. 

G/S  1 6.82  in. 

G/S  2 6.  96  in. 
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